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BEST AVAILABLE COPY 




SPECIFICATIONS /607 * 

1 . Name of invention 

Attitude controlling method in coaxial two wheeler 

2. Claim 

[ 1 ] With a coaxial two wheeler comprised of a pair of 
wheels (2, 3), axle (1) between the wheels (2, 3), car body (4) 
rotationally supported on the axle (1) , wheel driving motor (7) 
mounted on the car body (4), microcomputer (16) for controlling 
the wheel driving motor (7), and angle detection device (11) for 
detecting the angle of the car body (4); 

the slanted angle (01) of the car (4) detected by the angle 
detection device (11) is sampled at short intervals (At); the 
sampled angle (01) of the car body (4) used as conditional 
variables and feedback gain (K) used as a coefficient are 
respectively entered into the computer (16) to substitute and 
compute the pre-established input computation equation (60) to 
generate a control torque (T3) for the motor (7); and the 
computer (16) controls the motor (7) according to the computed 
control torque (T3). 

[2] With a coaxial two wheeler comprised of a pair of 
wheels (2, 3), axle (1) between the wheels (2, 3), car body (4) 
rotationally supported on the axle (1), wheel driving motor (7) 
mounted on the car body (4), microcomputer (16) for controlling 
the wheel driving motor (7), primary angle detection device (11) 

* Numbers in the margin indicate pagination in the foreign text. 
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for detecting the angle of the car body (4), and secondary 
detection device (14) for detecting the rotational angle of the 
wheels (2, 3); 

the slanted angle (01) of the car (4) detected by the 
primary angle detection device (11) and rotational angle (63) of 
the wheels (2, 3) detected by the secondary detection device (11) 
are sampled at short intervals (At); the sampled angles (01) of 
the car body (4), sampled rotational angle (03) of the wheels (2, 

3), and angle speed (01, 03) computed based on those samp-led 
angles (01, 03) are entered into the computer (16) as conditional 

variables. Those values (01, 03, 01, 03) and feedback gains 
-» -» 

(Kl, KT) used as coefficients are substituted in the pre-stored 

control equations (10, 49, 51) stored in the computer (16); a 

control torque (T3) of the motor (7) is computed based on this /608 

computation; the motor (7) is operated by the computer (16) 

according to the computed torque (T3). 

[3] With a coaxial two wheeler comprised of a pair of 
wheels (2, 3), axle (1) between the wheels (2, 3), car body (4) 
rotationally supported on the axle (1), attitude control arm (6) 
slantingly clamped to the car body (4), wheel driving motor (7) 
mounted on the car body (4), attitude control arm driving motor 
(9) mounted on the car body (4), microcomputer (16) for 
controlling the wheel driving motor (7) and. attitude control arm 
driving motor (9), primary angle detection device (11) for 
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detecting the angle of the car body (4), secondary detection 
device (14) for detecting the rotational angle of the wheels (2 , 
3), and third detection device (15) for detecting the attitude 
control arm angle (6); 

the slanted angle (61) of the car (4) detected by the 
primary angle detection device (11), rotational angle (03) of the 
wheels (2, 3) detected by the secondary detection device (11), 
and attitude control arm angle (62) detected by the third 
detection device (15) are sampled at short intervals (At); using 
the sampled angles (01)' of the car body (4), rotational angle 
(03) of the wheels (2, 3), slanted angle (02) of the attitude 
control arm (6), and angle speed (01, 03, 02) computed based on 
the sampled angles (01, 03) as conditional variables and feedback 
gains (Ki, Kt) as coefficients, those values (01, 03, 02, 01, 03, 
02) are substituted in the pre-entered control equations (20, 49, 
51) stored in the computer (16); control torques (T3, T2 ) for the 
wheel driving motor (7) and attitude control arm driving motor 
(9) are computed based on the values obtained from the equation 
computation, and the computer (16) controls the wheel driving 
motor (7) and attitude control arm driving motor (9) according to 
the control torques (T3, T2). 
3. Detailed explanation of this invention 
[Industrial field] 

This invention pertains to an attitude controlling method 
for a coaxial two wheeler. 




[Conventional technology] 

A coaxial two wheel vehicle having wheels at both ends of 
one axis can be made into more compact flat shape type device 
compared with a four wheel or three wheel car. However, it is 
not practically usable unless its attitude is stably controlled. 
Therefore, an attitude control method was developed, in which a 
control arm is hung from the car body supported on the drive 
shaft of the coaxial two wheeler to control the angle of the 
control arm against a car body. 
[Problems to be solved by this invention] 

However, with the attitude control method which drives a 
control arm in the direction opposite to the slanting direction 
of a car, not only does this method create a control delay caused 
by non-linear elements (e.g., non-sensitive area of the control 
arm driving motor, back-rash of the gear-transmission system, 
etc.), but also it is unable to eliminate the noise in the angle 
computation performed based on the short-interval sampling of car 
angles and control arm angles. Therefore, stable controlling on 
the car attitude using an actual feedback control is difficult. 
[Constituents of this invention] 
[Method to solve the problems] 

With the first claim of this invention, the tilt angle of a 
car detected by the first angle detection device is sampled at 
short intervals, and the control torque for the wheel driving 
motor is computed based on the sampled data, and the feedback /609 
control is performed on the wheel driving motor based on the 
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computed control torque. As a result, this makes the wheels to 
move in the slant direction of the car to provide a recovery 
moment to the car. With this method, despite the fact that the 
control torque is immediately computed from the sampled slant 
angles of the car, sufficient stable attitude control can be 
provided. 

With the second claim of this invention, the tilt angle of a 
car detected by the first angle detection device and wheel 
rotation angle detected by the secondary angle detector are 
sampled at short intervals, and the angle speed of those angles 
are computed. The control torque for the wheel driving motor is 
computed based on the sampled values and computed angle speeds to 
control the wheel driving motor using a feedback technique. As a 
result, the wheels can be operated in the slant direction of the 
car body to provide a recovery moment to the car. With this 
method, the use of angle speed computed from the sampled angles 
of the wheel and car body can further improve the attitude /610 
stability control. 

In the third claim of this invention, the tilt angle of a 
car detected by the first angle detector, wheel rotation angle 
detected by the secondary angle detector, and slant angle of the 
attitude control arm detected by the third angle detector are 
sampled at short intervals. Control torques for the wheel 
driving motor and attitude control arm driving motor are computed 
from the sampled values to control the wheel driving motor and 
attitude control arm driving motor using a feedback technique. 
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As a result, while the wheels are operated in the slant direction 
of the car body, the control arm rotates in the direction 
opposite to the slant direction of the car to provide a recovery 
moment to the car. With this method, the use of an attitude 
control arm can provide more improved control over the wheels and 
car body than the first and second methods. 
[Operational example ] 

The following explains the operational examples of this 
invention while referring to the figures. 

As shown in Figure 1, a pair of wheels (2, 3) are attached 
to both ends of an axle (1) and, in turn, a square frame-form car 
body (4) is tiltably supported on this axle (1). A fulcrum shaft 
(5) is rotating installed and supported in an upper part of the 
car body (4), and in the center of the fulcrum shaft (5), an 
attitude control arm (6) with a weight (6a) is slantably clamped 
thereto. A reciprocally rotatable car driving motor (7) is 
mounted on the car (4) immediately below the weight (6a). As 
shown in Figure 2, a deceleration gear array (8) is positioned 
between the driving shaft (7a) of the motor (7) and car axle (1), 
allowing the decelerated rotation of the wheel driving motor (7) 
to be transmitted to the wheel axle (1) to make the wheels (2, 3) 
reciprocatingly rotate. A reciprocatingly rotatable attitude 
control arm driving motor (9) is mounted on the car (4) 
immediately above the shaft (5), an array of decelerating gears 
(10) are placed between the driving shaft (9a) of the motor (9) 
and shaft (5). With this structure, the decelerated rotation of 
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the attitude control arm driving motor can be transmitted to the 
shaft (5) to guide the attitude control arm (6) in the 
backward/ forward directions. 

The primary rotary encoder (11) is attached to the side 
surface of the car (4), and the rotational shaft (11a) is placed 
on the extension line of the axil (1). A pair of contact pieces 
(12, 13) are diagonally mounted onto the rotational shaft (11a), 
and the tips of the pieces (12, 13) contact and engage onto the 
floor. With this structure, the tilt angle of the car (4) to a 
vertical line can be detected. The secondary rotary encoder (14) 
is mounted onto the wheel driving motor (7), and the third rotary 
encoder (15) is mounted onto the attitude control arm driving 
motor (9), detecting each turning angle of these wheels (2, 3) 
and a tilt angle to a vertical line. 

A microcomputer (16) for controlling the system is mounted 
at the lower area of the car body (4), and the detection signals 
transmitted from each rotary encoders (11, 14, 15) are inputted 
into the microcomputer (16) which computes the control torgues 
for the wheel driving motor (7) and attitude control arm driving 
motor (9) on the basis of these inputted signals, thus the motors 
(7, 9) can be controlled with the force eguivalent to the control 
torques. [1] Figure 3 shows the dynamic model of the simplified 
system exhibited in Figure 2. A wheel area (17) composed of a 
wheel part (19) consisting of an axle (1) and wheels (2, 3), 
motors (7, 9) supported by the car body (4), control computer 
(15), and car body (4), is expressed as a "material point" at a 
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distance (£1) from the axle (1); and an arm area (18) composed of 
an attitude control arm (6) and weight (6a) is expressed as a 
"material point" at a distance (£2) from the shaft (5). The 
wheel part (19) is considered to have the center of gravity on 
the axle (1), and the distance between the axle (1) and shaft (5) 
is £ . 

The motion energy (Ko) of the dynamic model is a total of 
motion energy (Kl) of the car body (17), motion energy (K2) of 
the arm part (18), and motion energy of the wheel part (19). The 
following describes the expression of each motion energy (Kl, K2, 
K3) using the coordinate elements of the x-y coordinate shown in 
Figure 2 where the mass of the car body (17) is expressed as ml, 
the mass of the arm part (18) is expressed as m2, and the mass of 
wheel part (19) is expressed as m3: 

Kl = (1/2) ml(xl 2 + yl 2 ) 
K2 = (1/2) m2(x2 2 + y2 2 ) 
K3 = (1/2) m 3 x3 2 

where x and y designate the first level of differential of x and y. /611 

When the tilts of the car part (17) and arm part (18) to a 
diagonal line are expressed as 01 and 02; the rotational angle of 
the wheel part (19) is expressed as 03; the semi-diameter of the 
wheels (2, 3) is expressed as r; inertia moment in the direction 
of axle (1) of the car part (17) is expresses as Jl; and inertia 
moment in the direction of shaft (5) of the arm part (18) is 
expresses as J2; the following expressions are formed: 
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K t - (1/2) m i 

x (( r*i 3 »iit -cos 0 ,)*♦ | t * 0 |* • sin 1 6 , } 

- (1/2) m i 

x ( r 1 8 3 J , 2 r I i 0 3 , -cos <? 1 1 I t * ^ 1 1 ) 

- (1/2) mi r 1 8 3 * 
+ mi r 2 £, 0 3 0, -cos 0|Ml/2) J| 0 1 1 

K2 - (1/2) m 2 
x (( r 8 3 i £ 8 t -cos tfi <j 0 2 -cos 8 2 ) * 
( - £ | • sin 0 i 2 0 2 . sin ,)) 

- (1/2) m 2 r 1 6*3* 

+ m 2 r^ 3 (/tfi • coj *i-£j «2 -cos 0 2 ) 
~ m 2 t i 2 & i & 2 • cos ( 9 i- 0 2 ) 

+ (1/2) m 2 I ' * iS(l/2) J 2 8 2 t 
K3 - (1/2) m j r » 0 3 » 

where 91, 02, and 03 are the first-level differential of 01, 02, 
and 03 . 

As for the motion energy (K 3 ) of the wheel part (19), 
rotational energy based on the inertia moment in the direction of 
the center of the gravity exists. When the inertia moment is 
designated as J 3 , the movement energy (K 3 ) can be expressed as 
below: 

K3 = (1/2) m3 r 2 03 2 + (1/2)J3 03 2 

Therefore, the motion energy (K) can be expressed as below: 

K o - K i ♦ K 7 + K a 
- (1/2) • ( J m 2 t ') 8 »»» (1/2) J 2 * 2* 
+ (1/2) *( m, r J tm 2 r * t m 3 r*» J 3 ) i ,» 

- m 2 t j t 8 , 8 7 • cos ( 8 1 - e 2 ) 

+ (m 2 <|»m 2 t ) r 8 x 0 j cos 6 , 

- m 2 I 2 r 8 2 *3 •' cos 8 2 

■ - ■ ( I ) 

The positional energy (P) can be expressed as below when the 
gravity acceleration ratio is expressed as g: 
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P » ( m i i t ♦ m 2 £ ) g - cos 5, 
- m 2 < 2 ff • cos 0 2 

• • • ( 2 ) 

The energy loss (D) can be expressed as below when the 
viscous friction coefficients of car-part motion (17), arm-part 
motion (18), and wheel-part rotation (19) are respectively 
expressed as f v f 2 , and f 3 : 

D » (1/2) f 2 * 2 1 + d/2) f 3 * 3" 

• • - ( 3 ) 

Since the torques (T2, T3) respectively exist for the arm 
part (18) and wheel part (19), a Lagrange motion equation 
described below can be established: 

d/dt*a/<Jtfl*K-<5/0*l*K 

(l - 1 . 2 . 3 ) 
- ■ ( 4 ) 

where d/dt is an ordinary differential equation, and 6/60i and 
6/60i are a partial differential equation. 

By substituting equations (1), (2), and (3) in Equation (4), 
and once Equation (4) is organized, the following equations can 
be established: 

At 0 i ♦ B t cos ( 0 , - 6 7 ) 6 1 s 
+ Ci cos 0, i'jtDi • sin ( 0,-0 2 ) 0 2* 
+ E| -sin ^ i • r t 0 t - 0 

- • • ( 5 ) 

Bi COS ( * 1 - 0 2 ) 0 | » B 2 0 2 
+ C 2 * cos 0 2 0 3»D 2 sin ( 0|-0 2 )0i* 
+ E 2 sin 0 f 2 02 -Tj 

- • - ( 6 ) 
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where 

iLj = + it^ Z 2 

B n = -iru, Z 2 Z 

= (ia, £ + nig Z)r 

D 1 = -nig jC 2 j0 

E 1 = - (n^ + m 2 Z )g 

B 2 = J 2 

C 2 = -n^ Z 2 r 

D 2 = -n^ Z 2 Z 

C 3 = . (in, + m 2 + m 3 )r 2 J 3 
D 3 = -(ia, Z + m 2 Z)r 
E 3 = ~m 2 Z z r 

[2] Next, the following explains the design of the linear type 
control system for the dynamic model shown in Figure 3: 
The following approximation is performed under the 
assumption that the dynamic model shown in Figure 3 is 
consistently controlled so that 61 and 02 never become 
excessively high values. 

J * n 9 I * d 1 tio « 2 w } 

cos 9 1 - cos 9 2 - cos ( 6 ,. 8 2 ) * { 

Jin 9 2 9 2 2 -sin ( 9 i . 0 z ) 8 ,« 

-sin ( 9 t . 0 2 ) e 2 * - 0 
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Based on this approximation, the equations (5), .(6)/ ar *d (7 
can be made into the following linear equations: 



0t ■ An« i#Au» j*A||5 iiAm* jiBiiTj 

• • ■ ( 8 ) 

it - Ail* itAji* 2 *Ajji M Ai4*i* BiiTa 

... ( 9 ) 

where 

A n - E i B j/A Aii--Ei Bi/A 

Ai3-fiB 2 /A A| 4 --r 2 Bi/A 
A j I ~ - E i B,/A A 22 - - E 2 Ai/A 
Ais- - f , B |/A A 2 4- - f 2 Ai/A 
B »«- Bi/A B 2 i- - At/A 

A " B i 1 - A 1 B 2 



Assuming that xl =01, x2 =02, x3 =01, and x4 = 02, the 
following vector and matrix display are performed: 



X " x 2 A - 0 0 0 1 

** A11A12A13A14 
x * 1 A21 A 22 A23A24 



9 1 r i 0 0 0 1 

inn' 



0 1 c~- 
B 11 
B 2 1 



u = T 2 



The following conditional equation and output equation can 
be obtained using the vector and matrix display described above: 



X ( t ) -X - X (t) *B-u <t> 

■••(10) 

Y ( t ) -*C • X ( t ) 

• • • (11) 



By Laplace-transforming the conditional equation (10), the 
following equation can be established: 



t u. -st 

X ( t ) c d t 
' o 

f Z -st , 1 -si 

- X(t)e d t + s I X ( i ) e dt 
' o o 

- s X ( s ) -X - X (s) + I • u { s ) 

Therefore, 1 

X ( s ) - ( s T-X) * B • u ( s ) 

■ • • (12) 

where I is a unit matrix, (si - A)~ 1 B is a reversed matrix of (si 
-A). 

By substituting the equation (12) in the output equation 
(11), the following equation can be established: 7613 

Y(s)-C~(sT-A~)B u(s) 

. • • (13) 

-> 

The function expressed as Y (s)/u(s) is called "transfer 
function" based on the control theory, and is expressed as the 
following polynomial equation: 

Y ( i ) / u ( s ) 
- k ( s - « 0 ( a - « a) ( * - * a) < » " * «> 

+ (J--It)( »-^3)( 

- • • (14) 

(s - Al), (s - A2), (s - A3), and (s - A4) in Equation (14) 
are characteristic polynomial equations, k is a gain, and Ai (i = 
1, 2, 3, 4) is called "pole". When the pole Ai exists on the 
right half of the complex variable plane, that is, the real 
values are positive, the system expressed with equations (10) and 
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(11) is known to be unstable; whereas when the pole Ai exists on 
the left half of the complex variable plane, that is, the real 
values are negative, the system expressed with equations (10) and 
(11) is known to be stable. For example, when the system 
parameters of the dynamic model shown in Figure 3 are the 
following values, the pole is confirmed to exist on the right 
half on the complex plane. 

ia, = 2.59 Kg £ = 0.318 m 

nig = 1.52 Kg i = 0.17 in 

m 3 = 0.36 Kg ^ = 0.20 in 

J 1 = 0.108 Kgm 

J 2 = 7.01 x 10 Kgm 

J 3 = 5.00 x 10 Kgm 

r = 4.5 x 10" 2 m 

^ = 6.0 x 10 Nms 
f 2 = 4.40 x 10 Kgm 
f 3 = 4.40 x 10 Kgm 

In order to stabilize the dynamic model shown in Figure 3, 

the control input u(t) must be added to transfer the unstable 

condition to the stable condition. To achieve this, the 

existence of u(t) must be guaranteed. If u(t) exists, the 

dynamic model expressed as Expressions (10) and (11) is 

controllable. It is a known factor based on the control theory 

-» -> 

that the necessary/sufficient condition is that the ranks [B, AB, 
Fb, ^B] of a matrix called controllable matrix match the 
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conditional variables (xl, x2, x3, x4). The system expressed 
with equations (10) and (11) is controllable as it satisfies this 
necessary/sufficient condition. 

To transform the system expressed as equations (10) and (11) 
to the system to which conditional feedback control is added, the 
following equation is used: 

«(D-K-X(t)+cV(t) 

■ • • (15) 

-» 

where K is a vector called conditional feedback gain, and c is a 
scaler . 

By substituting the equation (15) in equations (10) and 
(11), the following conditional equation and output equation can 
be obtained: 

X ( t ) - (A + BK) X (t) 
+ c B • V ( t ) 

• • • (16) 

Y C t ) -~c • X ( t ) 

• • • (17) 

By Laplus-transforming the conditional equation (16), the 
following equation can be made: 

X < i ) - (iT-I-BK) c B • V { 3 ) 
•■•(18) 

By substituting equation (18) to the output equation (17), 
the following equation can be made: 

Y ( i ) -"CdT-A-BK) c B • V ( 3 ) 
■ • • (19) 



c 

Therefore, inherent polynomial of the matrix (A + BK) 
matches 

the characteristic polynomial, thereby the inherent value 
of matrix (A + BK) agree with the pole. 

It is an established factor based on the control theory that /614 
equations (16) and (17) must be controllable for matching the 
characteristic polynomial expressed as equation (19) with Monique 
polynomial [i.e., necessary/sufficient condition for 
optionally establishing the inherent value of matrix (A + BK) ] . 
Therefore, by appropriately establishing the conditional feedback 

-» 

gain K, the system can be stabilized by positioning a pole at an 
appropriate location in the right half of the complex plane. 

When c = 0, the conditional equation (16) designating the 
controllable system shown in Figure 3 is expressed as below: 

X (t) - (X+BK) X ( t ) 

that is, 

u ( l ) - K ■ X < t ) • • • (20) 

The feedback gain K in equation (20) can be obtained by the 
following method: 

(1) Select preferable four different and symmetrical poles Si (i 
= 1/ 2, 3, 4). 

(2) Compute the following equation after selecting the real 
number (gi) when all the poles (Si) are real numbers. 
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Fi = (Si • I - A) Bgi 
(3) When all the poles (Si) are complex numbers a + flj (j is an 
imaginary number), Si + 1 is made into a conjugate root al - 
Bi j, and appropriate real numbers gi, gi+ 1 are selected, 
then, the following eguations are computed: 

fl - ( (or)T-X) B g It 0 ! BgUt ) 
+ ( ( oc I T- A~) *♦ 0 l*T) 

f hi - { (oriT-X) B glti^l B g I ) 

+ { ( of i T-T> 2 ♦ $ i*T) 



With those computations, the feedback gain K becomes the 
following: 

f r t i 

K - (g1.B2.g3.g4 J f 2 

f 3 
1 f 4 1 

-Cgl.g2.g3.g4) C f t . f 2. f 3. f 4 1 

• • • (21) 

[3] In special cases of the dynamic model shown in Figure 3, 
that is, the coaxial two wheelers from which the attitude control 
arm (6) is removed (shown in Figure 4 and its dynamic model is 
shown in Figure 5), The Lagrange motion equation of the dynamic 
model shown in Figure 5 can be obtained by the following equation 
using equations (5), (6), and (7): 

A i 0 C 1 - cos 8 t d 3 * E 1 -sin 0 x 
♦ t i 0 i - o 

■••(22) 

C I * C O S 6 i * i ♦ C 3 3 

+ D 3 sin 5, i,* f f 3 tf 3 „ Tj 
• ■ (23) 
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The following approximation is performed with the equations 
( 8 ) and ( 9 ) : 

With this approximation , equations (22) and (23) are linear- 
processed as described below: 

6 i - Au*i»Aui|*Ai4*i*BiiTj 
. • • (24) 

8 3 - A 21 e t » Aja 5 if A i* h j4 B 21 T 3 
--■(25) 

where 

An- Ei C 3 /A An - - f i Ca/A 

A 14- f 3 C i/ A An-Ei C|/A 

Au-f, Ci/A A 24 - - f a A|/A 

B | | - " Ci/A B 2 i - A i/A 

Ei - - mi g / 1 C 1 - m 1 < 1 r 

Cj - (mpinj) r 1 + J ] At - J i 
A - A ! Cj - C3 2 

Assuming that xl = 01, x2 = 03, x3 = 01, and x4 = 63, the 
following vector and matrix display are performed: 

— 1 x 1 1 _ t 0 0 1 0 * 

X - I x 2 I A-|0 0 0 1 

M A11A12A13A14 

1 x 4 ' 1 An A22 A 2 3 A 2 4 



0 1 r 1 0 0 0 

B 11 
B 2 r 



0 I ~C - [ 

B 11] 1 0 I 0 0 



u = T3 /615 
Using those vector and matrix display, the following 
conditional equation and output equation can be obtained: 
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X(t)-A-X(l)+B*tt(l) 
•(26) 

Y ( t ) - C" ■ X ( i ) 

• • • (27) 

Parts (a) and (b) in Figure 6 show the computer simulation 

-» 

results for the equations (26) and (27) using the feedback gain K 

expressed as expression (21) established according to the methods 
(1) ~ (3) described above. With those simulations, the numeric 
value differentiation method described later is used. In Part 
(a) in Figure 6, curve (CI) shows the variation of the tilted 
angle (61) of the car body part (17), and curve (C2) shows the 
variation of the rotational angle (03) of the wheel part (19). 
In Part (b) in Figure 6, curve (C3) shows the variation of the 
tilted angle speed (01) of the car body part (17); curve (C4) 
shows the variation of the rotational angle speed (03) of the 
wheel part (19). The parameters and conditions are established 
as following: 
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Note that X (0>= [0.0, tt/10.0] expresses the initial 
condition; disturbance n/10 (rad/s) is provided to the car part 
(17); and 5 ms is used as simulation interval At. The simulation 
results shown in Parts (a) and (b) in Figure 6 indicate that the 
car part (17) can be returned to the balanced position (01 = 0) 
shown in Figure 4 and maintained thereafter. Also, with this 
balancing control, although the further the pole positioned to 
the right half of the complex plane, the greater the control 
input becomes, the control system can be quickly stabilized. 
[4] With the dynamic models expressed as equations (22) and 
(23), assuming that the attitude control arm (6) is removed and 
the wheels (1) are almost stand still to ignore the motion 
related to the wheels (1) and approximation of sin 01 * 01 is 
performed, the following motion equation can be established: 

Ai 0 it f i * t* E 1 -stnfli 
= J l 8 it t i 0i-mi g/0i - 0 

• • ■ (59) 

Feedback of u(t) = T3 is provided to the system expressed a& 
Equation (59). 

J i 0 i ♦ I i 6 \ - m\ % t 6 | - u (t) 

- - e k e i c t > 

• • • (60) 



Note that this is a constant determined by the torque 
function of the wheel driving motor (7), deceleration ratio of 
driving system, etc. 
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Considering the overall control delay time aT [e.g., 
detection of tilt angle 01, computation time by the computer 
(16), command input time for wheel driving motor (7), transfer 
delay by the deceleration gears (8), etc.], the equation (60) 
becomes the following: 

J t £ i ♦ f i 8 i - ra i g t 6 i 

■-«K»i (t-AT) 
• • '(5 1) 

By Taylor-developing 01 (t - aT) while ignoring the two 
dimensional small values, the following equation is made: 

6 , ( t -AT) * d t < t ) - A T 6 i( l j /616 

As a result, equation (61) becomes the following: 
j , *, t (Ii-kut) i . 

+ (K{-mt g£) * I - 0 
... (62) 

The following characteristic equation can be obtained by 
Laplus-transf orming the equation (62): 

Ji * *♦ ( f i- K € A T ) s 

+ (K { - n»i | I ) - 0 

•••(63) 

The stabilization condition for the control system expressed 
as equation (63) can be obtained as below based on the Routh- 
Hurwitz stabilization conditions: 

nug*/f<K<eAT 

• • ■ (64) 
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With equation (64), although the lower limit of the feedback 
gain K is fixed, the upper limit varies depending on the delay 
time aT. That is, smaller the aT is, the greater the selection 
range of feedback gain becomes, and the selectability of 
appropriate feedback gain K becomes higher. 

The computer simulation result performed on equation (62) 
using feedback gain K selected from equation (64) exhibits the 
attitude stability control equal to the stable attitude stability 
provided to the linear control system described above (refer to 
[3]). For example, the attitude stability is excellent when aT = 
3ms and K = 10.0. 

Also, the experiment results on the control system expressed 
as equation (62) confirm that no negative results occur with the 
use of sampling time At when the delay time is aT. 
[5] The feedback controls described above are all based on the 
polar position method. The following explains the feedback 
control using an optimal control method: 

The control system shown in Figure 3 is expressed as 
equation ( 10 ) : 

X (t) - A~ • X (t) + B- u (t) 



The control input u(t) is then determined from equation (10) 
to minimize the evaluation function J(u) : 



(10) 




T 

+ u(l)RU(t))di 



(28) 



where Q and R are respectively 4x4 and lxl absolute symmetric 
matrixes. By minimizing J(u) as much as possible, the initial 
condition can be maximally stabilized with a small control 
energy. The optimal input uo(t) which minimizes J(u) can be 

obtained as equation (31) from the solution P of the Riccati 
equation (29 ) : 

T T 7 + TX -Tb R 1 b'T+q" - o 

•••(29) 

K R ' B? •••(30) 

u 0 ( t ) - K • X ( t ) ... (3l ) 



Therefore, by obtaining K, the optimum input uo(t) can be 
produced in the form of feedback. 

The feedback gain K for the conditional equation (16) 
expressing the controllable system shown in Figure (3) can be 
obtained in the following manner: 

(1) Appropriate diagonal matrix Q and R are selected. 

(2) K (i = 0) is selected to stabilize A + BK (K computed in the 
polar positioning method is used for simplification). 

(3) Assuming Ai = A + BKi, 4x4 semi-absolute symmetric matrix 
solution P of the Lyapnov equation (32) is obtained. 

Tl T + Tl At + RK-0 

• • (32) 
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(4) Return to (3) after setting Ki+1 = -R~ 1 B T Pi and~i = i + 1 . 

By repeating those processes, the convergence value K of Ki 
can be obtained. 

The weight Q and R in the evaluation function (28) strongly 
affect the control capability of the control system, and 
according to the results of experiments conducted by the 
developer of this invention, smaller the ratio of R to Q is, the 
faster the controlling speed becomes. Also, with the dynamic 
mode shown in Figure 3, since the stability of the car part (17) 
is more important than stability of the arm part (18), Q must be 
delibe-rately selected so that greater weight can be placed to 
the diagonal element of Q corresponding to XI. For this reason, 
the optimum control method is better than, the polar positioning 
method. As described with the simulation result obtained with 
the polar positioning method, The simulation result obtained from 
the optimum control method returns the car part (17) to the 
balanced position (01 = 0) (see a broken line in Figure 4), and 
the balanced condition can be maintained thereafter. 
[6] The following explains the design of a non-linear type 
control system of the dynamic model shown in Figure 3: 

With the linear type control system, an assumption was made 
that the dynamic model shown in Figure 3 was constantly 
controlled to prevent excessive 61 and 02. With the non-linear 
type control system, the following assumptions are made: 
(1) The slant angle 01 of the car part (17) is controlled so 
that it does not become excessively high. 
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(2) The rotational acceleration speed 03 of the wheel part (19) 
is extremely narrow. 

(3) When the slant angle 01 of the car part (17) is wide, the 
angle speed 01 is narrow. 

(4) When the slant angle 02 of the arm part (18) is wide, the 
angle speed 01 is narrow. 

(5) When the angle speed 01 of the car part (17) is wide, the 
slant angle 01 is narrow. 

(6) When the angle speed 02 of the arm part (18) is wide, the 
slant angle 02 is narrow. 

With those assumptions (1 - 6), equations (5), (6), and (7) 
expressing the dynamic model shown in Figure 3 can be simplified 
as described below: 

Ai 0|»B| -cos Or 0j»Ei tf|#fi 6\ - 0 

• • • (33) 

Bi -cos 0 t *i»B2 02*Ej • sin 0 t 
+ f 2 0 2 - T 2 

• • • (34) 

Ci dud -cos 0 r 0r*Ci 0 3 * f 3 0 3 - T 3 
••■(35) 

. Ox , J a b J , - Ei 0 t - f 1 0 ! 1 
I o\ I " ' c 4 I * T 2 ' 1 
.-•(36) 

/a b 1 , Ai B 1 • c o a 0 2 J -I 
'< d ' ' B 1 • c o s * 2 62 ' 

• • • (37) 

• Tj 1 - T 2 - Ei -sin Or - I 2 0 2 

•••(38) 
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Then, the following vector and matrix display are performed: 





r 6 1 




r 0 , 




e i 


B - 


b 








0 




1 h t J 




c ' 



A - 



0 111 

a E i - a f i 0 0 

0 111 

b E i — b f 2 0 0 



As a result, equation (34) is rewritten into the following 
conditional equation (39): 



e (t) - A • 0 (t) +b • T 2 ' (i> 

•••(39) 



Since A in equation (39) is a coefficient varies with the 
passage of time (time-variation), at least one pole is positioned 
at the right half of the complex plane, where the system 
expressed as equation (39) is known to be unstable. In order to 
fix the time-varying pole in the unstable system at a position in 
the left half of the complex plane to stabilize the system, the 
following feedback control is performed to the system expressed 
as Expression (39): 

T 

Tj" - tki.ki.kj.k4 1 '['i.'t.'j.«4 1 
-K8 

■ . • (40) 

As a result, expression (39) can be made into the following 
equation : 

e<t>-<A>BK)e<t) 

••(41) 
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The characteristic polynomial of equation (41) (i.e., the 
inherent polynomial equation A + BK) becomes the following: 
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D ( A) 

-dct [A T- (A + BK ) 1 

- I iT-X- BK I 

- >l 4 + ( a ( b k j- c k 4 ) A 3 

+ ( a E i- b k i- c ka- f i k J A ) >l * 
- ( f i ka* E i k «) J / A - E i k 3 / A 
• • • (42) 

where a = ac - b 2 . 

Since the poles in the system expresses as expression (41) 
are all assumed to be fixed, the coefficient in the equation 
(42) must 

not change with the passage of time (non-time-varying). Since A 

-> -> 
and B are time-varying coefficients, K naturally becomes time- 

-> 

varying. To obtain this time-varying K (from here, expressed as 
-» 

Kt), one preferably stable linear non-time-varying model 
expressed as Expression (43) is selected. 

X (t) -X ' • X ( t ) —-(43) 

As the control system expressed as equation (41) always 
agrees with equation (43) [i.e., the poles in the system 
expressed as equation (41) can agree with and is fixed to the 
poles in the system expressed as expression (43)], the control 
system expressed as expression (41) can be stabilized. The poles 
for the equation (43) are expressed as below: 

-al + 131 j (a > 0, i = 1, 2, 3, 4) 

This characteristic polynomial is expressed as below: 
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D ( i ) - I *T- A" * I 

X 0 0 0 
0*00 
0 0^0 
0 0 0 * 

- a 0 , j 0 0 0 
0 - cc ik $ 2 J 0 0 
0 0 - <*3> 0 ai 0 
0 0 0 - * 4* 9 4J 

- i U d 3 i d 7 * *» d i J + d o 

... ( 4 4 ) 



By consistently coincide the coefficient of the 

characteristic polynomial in equation (42) with the 

characteristic polynomial of the equation (43) , each element (kl, 
-> 

k2, k3, k4) of Kt becomes the following: 

k t - 1 / b ♦ ( a E i - d 2 ♦ f i <* t / E i 

-c Ado/Ei-d'do/Ei 1 ) 
• • • (45) 
k 2 - 1/ b ■ < a I|.dj»c Ad|/Ei 
+ f i c Ado/Ei 1 ) 

-♦•(46) 

k 3 - -Ado/Ei 

■ • ■ (47) 

k t - (A/Ei)-( fi do/Ei-di ) 
•••(48) 

The poles in the equation (39) expressing the time-varying 

system produced by those computations consistently coincides with 

the poles (-al + fllj) of the stable linear time-varying model 

-> 

(41), and the feedback coefficient Kt is made into a time-varied 
quantity at each sampling time by equations (45) - (48). That 
is, this feedback becomes a kind of adaptation-control process, 
which is a model-matching process to the time-varying 
system parameters. 
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The control input Tl of the arm part (18) can-. be obtained 
from equation (38) by the following method: 

Tj - T 2 ' + E: -sin 0 2 ♦ f 2 d 2 
- Kt 6 + E 2 '5»n * 2» f 2 8 2 

• • (49) 



The control input T3 for the wheel part (19) can be obtained 
from equation (35): 

Tj.C, *|tC| 'COS 8 2 i?2tC3 *3»f3 *3 

= [c,.c 2 -cos ) + '3 ij»c 3 9l 

By substituting equations (37) and (40), the following 
equation is made: 

- [ c t , C 2 ■ cos 9 2 ) 
-Ei * |- ( i 6 i 



a b 



b c 



r -E|ff|-I|Vli , • 

I Kt • e J 



• ■• - (50) 

By entering a control input that can approximate the angle 
speed 03 of the wheels (2, 3) close to the target angle speed 03, 
the wheels (2,_ 3) can be transferred at equal speed. Therefore, 
the control input for the wheel part (19) is expressed as below: 



T 3» ( C i . C 2 * COS 0 2 1 



a b 



b c 



-Ei d i- i t $ » 

K • e 



+ f 3 0 a C 3 ( 6 3- 0 3 • ) 

• • - (51) 
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The time-varying feedback gain Kt is computed from the 
sampling condition quantity (01, 02, 03) of each time. Parts 
(a), (b), and (c) in Figure 7 shows the computer simulation 
result performed on equation (39) using those computed values. 
The curves (C5, C7, C9 ) shown in Parts (a), (b), and (c) in 
Figure 7 designate the variation of the slant angle 01 of the car 
part (19), and curves (C6, C8, CIO) designate the variation of 
the slant angle 02 of the arm part (18). The parameters and 
conditions for those simulations are the following: 
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Note that the X(0) = [0, tt/10, 0.0] expresses the initial 

condition, and a disturbance tt/IO (rad/s) is provided to the car 
part (17), where the duration of 10 ras is used as simulation 
interval At. 

The pole (Si) in Part (a) in Figure 7 is expressed as: 



Si = 



(-3±4j, - 



4±5j) (i = 1, 2, 3, 4) 



The pole (Si) in Part (b) in Figure 7 is expressed as: 



Si = 



(-5±4j, - 



6±5j) 
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.The pole (Si) in Part (c) in Figure 7 is expressed as: 
Si = (-7±4j, -8+5j) 

The simulation results shown in Parts (a), (b), and (c) in 
Figure 7 indicate that the car part (17) can be returned to the 
balanced position (01 = 0) (see a chained line shown in Figure 2) 
and maintained at the position thereafter. Also, with this 
balancing control, although the further the pole is positioned to 
the right half of the complex plane, the greater the control 
input becomes, the control system can be quickly stabilized, /620 

The computer simulation results of system designs described 
in [2] - [6] all exhibit the control effectiveness of the dynamic 
model shown in Figure 3 and Figure 4. The coaxial two wheeler 
shown in Figure 2 is formed based on this theoretical argument. 
The schematic of the control system for the coaxial two wheeler 
is shown in Figure 10, and a control program shown in Part (a), 
(b), or (c) in Figure 11 is stored in the computer (16). The 
control program shown in Parts (a) and (b) in Figure 11 
correspond to the linear type control system, whereas the control 
program shown in Part (c) in Figure 11 corresponds to the non- 
linear type control system. 

With the linear type control program shown in Part (a) of 
Figure 11, the computer (16) initializes various control circuits 
(e.g., counter circuit, output interface circuit, data latch 
circuit, etc.) when a control initiation signal is inputted, and 
processes detection signals outputted from the rotary encoder 
(11) via a counter circuit. The computer then computes the 
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control torque (T3) using the equation (60) consisting of the 
sampled value 01 and pre-established feedback gain K, and outputs 
the computed control torque (T3) to the drive part via a data 
clutch circuit and D/A converter shown in Figure 10. The data 
clutch circuit holds the output data transmitted from the 
computer (16) until the next sampling point is reached. 

With the linear type control program shown in Part (b) in 
Figure 11, the computer (16) initializes various control circuits 
(e.g., counter circuit, output interface circuit, data latch 
circuit, etc.) when a control initiation signal is inputted, and 
simultaneously processes the detection signals outputted from the 
rotary encoders (11, 14, 15) via a counter circuit. The computer 
then computes the differential values of the nth sample (01, 02, 
03) selected from the n groups (01, 02, 03) (i) (i = 1, 2...n) 
[i.e., angle speeds (01, 02, 03) (n)]. 

The computer then computes the control torque (T2 and T3) 
using the equation (20) comprised of sampling values (01, 02, 03) 
(n), computed values (01, 02, 03) (n), and pre-established 
feedback gain K. The computed control torques (T2 and T3) are 
output to the drive part via the data latch circuit and D/A 
converter shown in Figure 10. 

With the non-linear type control program shown in Part (c) 
in Figure 11, the computer (16) performs the same process as 
described in the linear type control program up till the angle 
speed computation process using numeric value differentiation, 
and then computes the time-varying feedback gain Kt. The 
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computer (16) computes the control torques T2 and T3 using 
equations (49, 51) comprised of sampling angles (61, 02, 83) (n), 
computed angle speed (61, 02, 03) (n), and pre-set time-varying 
feedback gain Kt, and transmits those torques (T2, T3) to the 

drive part via the data latch circuit and D/A converter shown in 
Figure 8. 

The following explains the numeric value differential method 
performed by the linear and non-linear control programs: 

Sampling time At [i.e., between the incorporation times of 
(01, 02, 03) (1) and (01, 02, 03) (t+1)] is sufficiently small, 
the angle 0i (i = 1, 2, 3) can be approximated by a two- 
dimensional polynomial described below: 

0i(t) = a + bt + ct 2 

... (52) 

By differentiating the equation (52), the angle speed 01(t) 
can be obtained in the following manner: 

0i(t) = b + 2ct 

. . . (53) 

The sum of squares of the difference of observation angles 
0i(l), 0i(2), 0i(3), 0i(4), and 0i(5) at the n pieces of sampling 

points (5 pieces in this example) 
tl = -At t2 = -At t3 = 0 
t4 = At t5 = 2At 

and 0i can be expressed as below: 



34 



o 
r 



& = [0i(l) - 0i(l)] + [6i(2) - 0i(2)] 
+ [0i(3) - 0i(3)] + [0i(4) - 0i(4)] 
+ [0i(5) - 0i(5) ] 

- I fl i (II - ( a + 2 b A. t + 4 c A t ") } 1 
+ ( 0 | (2) - <a + bAt+cAi*))» 
+ ( 8 i (3) - a ) * 

+ ( * i (4> - ( a + b A t + c A t ') > ' 

+ C ^ i <S) — (a + 2bAt+4cAt 1 ))* 

The coefficient a, b, and c of the equation (52) are 
determined in the following manner to minimize the value. of "e" 

a - U/70) - (-6 d t (0*24 8 \ (2) 

♦ 34 8 | (3U24 0 \ W). 6 9 <S>) 

b = (1/10 At) . <- 2 * i (I) - 6 \ (21 

♦ e j (4) ♦ 2 * i (sj) 

c - (1/14 A t * ) . (2 0 i tl) - 8 \ (2J 

• 2 * I (31 - * i (4) -2 0 i (5)) 

By substituting the a, b, and c determined in this manner 
the equation (53), the angle speed 0i(n) at each sampling point 
(n = 1, 2, 3, 4, 5) can be obtained. For example, the angle 
speed at the sampling point (t) can be expressed as below: 

i i (51 ■ b + 2 c • 2 A t 

- (1/10 At) • ( 26 0 i (1) 27 6 i (2) 

-40 8 \ (3)- 13 0 i M)»54 0 i (51 ) 
♦ • • (54) 

That is, the angle speed 0i at the sampling point (t) is 
computed while the angle ©i at the sampling point is smoothed 
based on the past five sampling angles 0i(l) -01(5), and the 
computed angle speed wi^h a smoothing effect precisely agrees 
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with the actual angle speed. 

The length of the sampling time interval At affects the 
stability of the control system. The smaller the At is, the 
better the results become, where the computation speed of the 
computer (16) between two sampling points determines the interval 
At. Therefore, the equations (45) ~ (48) for computing each 
element of the time-varying feedback gain Kt in the non-linear 
type system are rewritten in the following manner: 

k i - (1/ b ) • ( a B i - d t+ f I d i/ B i 
- c A do/Bt- f t'do/Ei 1 ) 

- (1/ b ) - ( a E i - d 2 - f i k 4/ A - c k 3 ) 

••(55) 
k 2 - (1/ b ) • ( a f i- d s* c A d |/E i 
+ f t c A do/Ei 1 ) 

- (1/ b ) • ( a f t- d j. c k 4 ) 

••(56) 

k 3 A d o/ E i 

♦••(57) 

k 4 - (A/E|) • ( f , do/Ei-d] ) 
■ - (l/Ei)'( d| o + fj k] J 
•••(58) 

Using equations (55) - (58) instead of equations (45) ~ 

(48), 

-» 

the time-varying feedback gain Kt computation time can be 
drastically reduced. 

Parts (a), (b), and (c) in Figure 8 show the feedback 
control results for the coaxial two wheeler shown in Figure 4 
where the wheeler was controlled using the system parameter 
values and initial disturbance identical to those used with the 
linear type control simulation performed on the dynamic model 
shown in Figure 5, and operated at conditional feedback gain K = 
[4.88, 1.31, 5.37, 6.00] obtained for each pole Si = (-3+ j, -4± 
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2j). The curve (Cll) shown in Part (a) in Figure 8 designates 
the variation of the slant angle 61 of the car part (17), and 
curve (C12) designates the variation of the rotational angle 03 
of the wheel (19). The curve (C13) shown in Part (b) in Figure 8 
designates the variation of the angle speed 01 of the car part 
(17), and curve (C14) designates the variation of the control 
torgue (T3). The results shown in Parts (a) - (c) in Figure 8 
indicate strong stability of the system during the movements of 
the wheels (2, 3) allowing the car part (17) to immediately 
regain the balanced position after experiencing a significant 
disturbance. 

Parts (a) and (b) in Figure 9 show the feedback control 
results when a coaxial two wheeler shown in Figure 2 operated at 
the target of 03 = 3 rad, where the system parameters were egual 
to those used in the non-linear type control simulation for a 
dynamic model shown in Figure 3. The curve line (C15) in Part (a) 7622 
in Figure 9 shows the variation of the tilt angle 01 of the car 
body (17), curve line (C16) in Parts (a) and (b) in Figure 9 
shows the variation of the rotational angle 03 of the wheel part 
(19), and curve line (C17) in Part (b) in Figure 9 shows the 
variation of the angle speed 01 of the car body (17). In Parts 
(a) and (b) in Figure 9, as the car part (17) is slightly 
fluctuating at the beginning, the control process for moving the 
wheels (2, 3) is performed in the direction preferred for 
stabilizing the attitude of the car body (i7), controlling the 
wheels at short intervals. The arm part (18) also slightly 

37 



fluctuates in the same phase as that of the car body part (17) 
after the wheels (2, 3) start moving, and the fluctuation 
movements of the arm part (18) and car part (17) become 
synchronous approximately 7 seconds after the initiation of the 
wheel movement, causing the arm part (18) to greatly fluctuate. 
At this point, the wheels (2, 3) rapidly get close to the target 
position. The stability of the car part (17) during the 
acceleration process suppress the fluctuation of the arm part 
(18). The wheels (2, 4) starts decelerating at approximately 9th 
second. The instability of the car part (17) caused by this 
deceleration is compensated by the intermediate fluctuation of 
the arm part (18). The coaxial two wheeler stops slightly before 
the target position 03 = 3 rad. This is probably because of the 
friction of the solid parts of the wheel driving system. 

The control effect of this device when the target position 
is set to 03 = 0 was identical to the condition after the car 
stopped moving [see Parts (a) and (b) in Figure 9]. 

When the fluctuation of the car part (17) during the 
motionless period of the wheeler was compared between Part (a) in 
Figure 8 and Part (a) in Figure 9, the attitude control method 
which simultaneously drives the wheel part (19) and arm part (18) 
could suppress the fluctuation of the car part (17) far better 
than the method which could only control the attitude. 
Therefore, the coaxial two wheeler equipped with an attitude 
control arm (6) shown in Figure 2 is proven to be more effective 
in controlling the attitude of the wheeler than the coaxial two 
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wheeler not equipped with an attitude control arm. 

The developers of this invention conducted an attitude 
control experiment using only the attitude control arm (6). As a 
result, neither the linear or non-linear type control system 
could stably maintain the car body (17) at a balanced position. 
This is probably due to the control delay caused by the non- 
linear type elements (e.g., non-sensitive area of the motor (9) 
for driving an attitude control arm or back-lash of the gear 
transmission system, etc.) and the noise occurring during the 
angle speed computation performed based on the short-interval 
sampling of the slant angle of the car and attitude control arm 
(6) utilized with feedback control, revealing the difficult 
adaptation of the conventional attitude control method using only 
an attitude control arm (6). 

This invention is not limited to the operational example 
described above, and any modifications [e.g., providing a 
steering capability to the coaxial two wheeler by connecting the 
wheels (2, 3) and wheel driving motor (7) via a differential 
movement gear system, a pair of wheels (2, 3) are independently 
driven to provide a steering function]. 
[Effectiveness of this invention] 

With the first claim of this invention, the tilt angle of a 
car detected by the first angle detection device is sampled at 
short intervals, and the control torque for the wheel driving 
motor is computed based on the sampled data, and the feedback 
control is performed on the wheel driving motor based on the 
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computed control torque. As a result, this makes the wheels to 
be shifted in the slant direction of the car body to provide an 
immediate positional recovery of the car body. 

With the second claim of this invention, the tilt angle of a 
car detected by the first angle detection device and wheel 
rotation angle detected by the secondary angle detector are 
sampled at short intervals; the sampled angles of the car body, 
rotational angle of the wheels, and angle speed computed based on 
those sampled angles are entered into the computer as conditional 
variables. Those values used as control variable and feedback 
gains used as coefficients are substituted in the pre-stored 
control equations stored in the computer to generate a control 
torque of the motor; the computer can control the motor by 
feedbacking the computed torque. As a result, when the car body 
is slanted, the wheels can be operated in the same slant 
direction of the car body to assure the car body to regain the 
original position. 

With the third claim of this invention, the tilt angle of a 
car detected by the first angle detector, wheel rotation angle 
detected by the secondary angle detector, and slant angle of the 
attitude control arm detected by the third angle detector are 
sampled at short intervals. The pre-stored control input 
computation equations are computed by substituting the sampled 
angles of the car body, rotational angle of the wheels, slanted 
angle of the attitude control arm, and angle speed computed based 
on the sampled angles as conditional variables and feedback gains 
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as coefficients. The control torque.s for the wheel driving motor 
and attitude control arm driving motor are computed based on the 
equation computation described above, and the computer controls 
the wheel driving motor and attitude control arm driving motor 
according to the control torques, allowing the motors to be 
feedback-controlled. As a result, when the car body is slanted, 
the wheels are operated in the same direction of the slanted car 
body, while the control arm is rotated in the direction opposite 
to the slant direction, assuring the car body to regain the 
original attitude. With this method, the use of an attitude 
control arm can more assuredly maintain the balanced position of 
the car body than the first and second methods. 
4. Simple explanation of the figures 

Figure 1 shows the diagram of the device used to explain the 
third claim of this invention, and Figure 2 shows the cross- 
sectional diagram of the device shown in Figure 1. Figure 3 
shows the diagram of a dynamic model of simplified system 
displayed in Figure 2. Figure 4 shows the diagram of the side 
surface of the device used in the operational example. Figure 5 
shows the diagram of the side surface of a dynamic model. Parts 
(a) and (b) in Figure 6 show the computer simulation results of 
the dynamic model in Figure 5, where Part (a) shows a graph 
exhibiting the variation of the slant angle 01 and rotations 
angle 03 of the wheel (19), and Part (b) shows a graph exhibiting 
the variation of the angle speed 0 and rotational angle speed 03. 
Parts (a), (b), and (c) in Figure 7 shows the computer simulation 
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results performed on the dynamic model in Figure 3, where part 

(a) , (b), and (c) respectively show the graph exhibiting the 
variation of slant angles 01 and 02 of car part 17 and arm 18 
based on the selected pole Si = (-3+4j, -4±4j), pole Si = (-5±4j, 
6±4j), and pole Si = (-7±4j, 8±4j). Parts (a), (b), and (c) in 
Figure 8 show the feedback control results of the coaxial two 
wheeler shown in Figure 4; where part (a) shows the graph 
exhibiting the variation of slant angles 01 and rotational angle 
03 of wheel (19); part (b) shows the graph exhibiting the 
variation of rotational angle 03 of wheel (19); and part (c) 
shows the graph exhibiting the variation of control torque. 
Parts (a) and (b) in Figure 9 show the feedback control results 
of the coaxial two wheeler shown in Figure 2; where part (a) 
shows the graph exhibiting the variation of slant angles 01 and 
rotational angle 03 of wheel (19); part (b) shows the graph 
exhibiting the variation of rotational angles 03 and slant angle 
02 of the arm part (18). Figure 10 shows the diagram of the 
coaxial two wheeler shown in Figures 1 and 2. Part (a) in Figure 
11 shows the flow chart of a linear type control program. Part 

(b) in Figure 11 shows the flow chart of another linear type 
control program. Part (c) in Figure 11 shows the flow chart of 
non-linear type control program. 

1... Axil; 2, 3... Wheel; 4... Car body; 6... Attitude 
control arm; 6a... Weight; 7... Wheel driving motor; 9... 
Attitude control arm driving motor; 11... First angle detection 
device rotary encoder; 14... Second angle detection device rotary 
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encoder; 15... Third angle detection device rotary encoder; 16... 
Computer; 17..'. Car body; 18... Arm; 19... Wheel part; 91, G2 , 
03... Angle; 01, 02, ©3... Angle speed; T2, T3... Control torque; 



K, K, Kt... Feedback gain 
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Figure 7 (a) 

Key: a) Angle (rad); b) Time (S). 
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Figure 7 (c) 

Key: a) Angle (rad); b) Time (S). 
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Figure 8 
To 3 



Key: A) Angle (rad); B) Time (S); C) Angle speed (rad/S); D) 
Torque (Nm) . x ' 9 ' 9 
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Figure 9 (a) 

Key: A) Angle (rad); B) Time (S) 
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15) 



16) 



Figure 10 
Key: 1) Drive part; 2) Motor 9; 
Motor 7; 6) Drive part; 7, 8, 9) 
12) Control computer 16; 13) D/A 
circuit; 15) Data latch circuit; 
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3) Encoder 15; 4) Encoder 14; 5) 
Counter circuit; 10) Encoder 11; 
converter; 14) Data latch 
16) D/A converter. 
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5) 



Figure 11 (a) 

Key: 1) Start; 2) Control circuit initialization; 3) 01 
incorporation; 4) Compute control torque T3 from equation (60) 
comprised of ei and established feedback gain K; 5) Computation 
control torque T3 output. 
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6) 



Figure 11 (b) 

Key: 1) Start; 2) Control circuit initialization; 3) (01, 02, 
and 03) (i) incorporation; 4) Compute (01, 02, and 03) (n) using 
numeric differentiation; 5) Compute control torques T2 and T3 
from equation (20) comprised of (01, 02, and 03) (n) and (01, 02, 
and 03)(n) and established feedback gain K; 6) Computation 
control torque T2 and T3 output 
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6) 



7) 



Figure 11 (c) 

Key: 1) Start; 2) Control circuit initialization; 3) (01, 02, 
and 03) (i) incorporation; 4) Compute (01, 02, and 03) (n) using 
numeric differentiation; 5) Time-varying feedback gain Kt 
computation; 6) Compute control torques T2 and T3 from equation 
(49 and 51) comprised of (01, 02, and 03) (n) and (01, 02, and 
03) (n) and established feedback gain K; 7) Computation control 
torque T2 and T3 output. 
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